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' We report results of systematic numerical analysis of collisions between two and 

Qh. three stable dissipative solitons in the two-dimensional (2D) complex Ginzburg- 
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Landau equation (CGLE) with the cubic-quintic (CQ) combination of gain and loss 
terms. The equation may be realized as a model of a laser cavity which includes the 
' spatial diffraction, together with the anomalous group-velocity dispersion (GVD) 
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and spectral filtering acting in the temporal direction. Collisions between solitons 
are possible due to the Galilean invariance along the spatial axis. Outcomes of the 
collisions are identified by varying the GVD coefficient, (3, and the collision "veloc- 



. ity" (actually, it is the spatial slope of the soliton's trajectory). At small velocities, 



two or three in-phase solitons merge into a single standing one. At larger velocities. 



I both in-phase soliton pairs and pairs of solitons with opposite signs suffer a transi- 

■ tion into a delocalized chaotic state. At still larger velocities, all collisions become 

O ■ 

quasi-elastic. A new outcome is revealed by collisions between slow solitons with 
^ ' opposite signs: they self-trap into persistent wobbling dipoles, which are found in 

two modifications - horizontal at smaller /?, and vertical if (3 is larger (the horizontal 
ones resemble "zigzag" bound states of two solitons known in the ID CGL equation 
of the CQ type). Collisions between solitons with a finite mismatch between their 
trajectories are studied too. 



I. INTRODUCTION 



Complex Ginzburg-Landau equations (CGLEs) constitute a vast class of models for the 
pattern-formation dynamics and spatiotemporal chaos in one- and multidimensional non- 
linear media combining dissipative and dispersive/diffractive properties In particular, 



stable localized pulses ( "dissipative solitons" Q]) can be supported by CGLEs that meet the 
obvious necessary condition of the stability of the zero background. This condition rules 
out the simplest cubic CGLE, whose one-d^ensional (ID) variant admits well-known ex- 
act analytical solutions for solitary pulses [Si]. The stability can be achieved in systems of 
linearly coupled equations, with one featuring linear gain and the other - linear loss [4]. In 
such a model, exact stable solutions for ID solitons are available 5|]. Another possibility is 
to use the CGLE with the cubic-quintic (CQ) combination of nonlinear terms. For the first 
time, the CGLE of the CQ type was introduced by Petviashvili and Sergeev 6j in the 2D 
form, with the intention to construct stable fully localized 2D states. In ID, stable dissipa- 
tive solitons of the CQ CGLE had been later studied in detail [7|, including the analysis of 
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two-soliton bound states 
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121 Then, stable fundamental solitons 



lOj and localized 

vortices (alias spiral solitons) [lO|, [l^ have been found in 2D and 3D 1^ models of the 
CQ-CGLE type, as well as in the Swift-Hohenberg equation with the CQ nonlinearity 16|. 
Such equations find their most significant physical realization as models of large-area laser 



cavities, where the CQ combination of the 



linear amplifiers and saturable absorber_s 
In most above-mentioned works 



OSS and gain is provided by the integration of 



6j, [l0|-[16|], localized pulses and vortices were obtained 
as solutions to isotropic 2D equations. On the other hand, the CGLE which governs the 
spatiotemporal evolution of light in the large-area laser cavity is anisotropic, as its includes 
"diffusion" (the spectral filtering) acting only along the temporal variable. The existence of 
stable fully localized pulse solutions in the latter case suggest a possibility of the experimental 
creation of "light bullets", i.e., spatiotemporal optical solitons, in the cavities. In other 
physical contexts (unrelated to optics), anisotropy of the 2D CGLE was introduced in a 
different form, through unequal diffusion coefficients in the two perpendicular directions 

Q. 

In Refs. [91], stable spatiotemporal dissipative solitons were found in the model of the 



laser-cavity type, based on the following normalized CGLE with the CQ nonlinearity: 

tUz + \uxx + \{(3-t)UTT = - \tU + (1 - Z7i) |[/pf/ + i-i2\U\^U] . 



(1) 



Here, Z and X are the propagation and transverse coordinates in the cavity, and T = t—Z/Vo 
is, as usual, the reduced time, with t the physical time and Vq the group velocity of the 
carrier wave. Term Uxx in Eq- ([I]) represents the transverse diffraction in the paraxial 
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approximation, the coefficients accounting for the above-mentioned spectral filtering, Kerr 
nonlinearity, and background linear loss are all scaled to be 1, while j3 > corresponds 
to the group- velocity dispersion (GVD). Usually, a necessary condition for the existence of 
temporal solitons is (3 > [l^ , which implies the anomalous type of the GVD (in the present 
model, spatiotemporal solitons also tend to be more stable at /3 > jol). Further, positive 
coefficients 71 and 72 in Eq. ([1]) account for the cubic gain and quintic loss, respectively, 
which are characteristic features of CQ models. The third-order GVD was also taken into 
regard in Refs. j^, but this term is not considered here, as it does not essentially affect the 
results reported below. Because it combines the diffiaction along X and effective diffusion 
along T, Eq. ([T]) is called the diffractive-diffusive CGLE 91]. 

Once 2D solitons are available, an issue of obvious interest is to explore collisions between 
them, provided that they are mobile, i.e., the equation is Galilean invariant. The 2D CGLE 
with no diffusion obviously satisfies this condition, allowing free motion of solitons or local- 



ized vortices in any direction. This property was used in Ref. [lOj] to study collisions between 
solitons in the isotropic CQ CGLE, as well as their motion in external potentials. It was 
concluded that collisions between fundamental solitons result in their quasi-elastic passage 
through each other (with a resultant increase of the relative velocity), or mutual destruction 
of the solitons, or their merger into a single 2D pulse. In the same model, collisions between 
vortices demonstrated a quasi-elastic rebound. 

The laser-cavity model based on Eq. ([T]) features the Galilean invariance along the X- 
direction, which means that a moving solution can be generated from a quiescent one by the 
application of the Galilean boost corresponding to arbitrary "velocity" P (in fact, P is the 
tilt in the {X, Z) plane): 

U (X, T, Z) exp [i [PX - P^Z/2)] U {X - PZ, T, Z) . (2) 

This possibility suggests to consider collisions of 2D solitons in this model too. In this work, 
we report results obtained by means of systematic simulations of collisions between two and 
three solitons in the framework of Eq. ([1]). In the former case, both head-on collisions and 
those with a finite offset {aiming distance) between trajectories of the two solitons will be 
studied. In either case, we consider collisions between in-phase and out-of-phase 2D solitons 
(the latter means that they have opposite signs). 

In Section II, we report the results for two-soliton collisions, and in Section III - for 
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interactions between three solitons. Outcomes of the colhsions between two in-phase sohtons 
include the quasi-elastic passage at large velocities, delocalization in the X-direction (merger 
into an expanding quasi-turbulent state) at intermediate velocities, and merger of slowly 
moving solitons into a single stable pulse. A major difference for collisions between out-of- 
phase solitons is that, at small velocities, they do not merge into a single pulse; instead, 
they may form a new localized object - a wobbling dipole, i.e., a robust bound state of two 
solitons with opposite signs, which feature persistent oscillations relative to each other in the 
spatial direction. Moreover, two different species of the wobbling dipoles are reported below, 
horizontal and vertical ones. In the latter case, the out-of-phase solitons, although they 
collide head-on, shift in opposite perpendicular directions (along the T-axis), and eventually 
form a dipole with a fixed vertical separation between them. Unlike the results for collisions 



between dissipative solitons in the 2D isotropic CGLE [10|, in the present model we have 
never observed complete destruction (decay) of colliding solitons. For collisions with a finite 
aiming distance AT, we identify a critical value of AT which separates interactions and 
the straightforward passage. Three-soliton configurations feature either merger into a single 
pulse, or the transition into a delocalized chaotic state. In terms of the optical cavities, the 
various outcomes of the collisions offer possibilities for the use in all-optical data-processing 
schemes. 



II. TWO-SOLITON COLLISIONS 

A. The numerical procedure 

Equation ([T]) was solved by means of the 2D split-step Fourier method with 256 x 256 
modes and periodic boundary conditions in X and T, for the fixed size of the integration 
domain in both directions, \X,T — 10| < 10. The stepsize for the advancement in Z was 
0.01. To generate the first stable 2D pulse boosted to velocity (tilt) P, cf. Eq. ([2]), an initial 
configuration was taken as 

Uo (X, T) = exp [- (X2 + t2) /4 + tPX] , (3) 

see the first panel in Fig. [2] below. The numerical integration of Eq. ([1]) led to quick self- 
trapping of the input pulse into a moving (tilted) dissipative soliton, which is an attractor 
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of the model. The profile of the established soliton can be seen in the first panels of Figs. [3] 
and [51 

Generic results for collisions between the solitons with velocities ±P can be adequately 
represented by fixing the cubic gain and quintic loss coefficients to be 71 = 2.5, 72 = 0.5, 
while varying P and GVD coefficient (3. To generate diagrams presented below in Figs. [H 
and m which display outcomes of the collisions, we changed P and P by small steps, the 
initial configuration for each simulation being a stable pulse produced by the simulations at 
the previous step. 



B. Head-one collisions between in-phase solitons 

Outcomes of collisions between two identical stable solitons, set by kicks ±P on the head- 
on collision course, are summarized in Fig. [H Stable solitons exist only for (3 > /5min ~ —0.5, 
which determines the left-hand edge of the diagram. 

The simplest outcome of the collision is the straightforward quasi-elastic passage of the 
solitons through each other. We do not illustrate it by a separate picture, as it seems quite 
obvious; as well as in Refs. 9|, the solitons keep the mutual symmetry after the quasi- 
elastic collision, and demonstrate some increase of "velocity" P (recall it is actually defined 
as the tilt of the soliton's trajectory in the (X, Z) plane). Moreover, running the simulations 
in the domain with periodic boundary conditions, we observed multiple quasi-elastic colli- 
sions between solitons. The solitons which emerge unscathed from the first collision survive 
indefinitely many repeated collisions as well. 

With the decrease of the collision velocity, the quasi-elastic passage is changed by the 
delocalization. This means that two in-phase solitons, interacting attractively, merge into a 
single pulse, which, however, fails to self-trap into a standing soliton. Instead, it gives rise 
to a quasi-chaotic ("turbulent") state, that remains localized in the temporal direction (T), 
but features indefinite expansion along X, see a typical example in Fig. [2l This outcome 
may be explained by the fact the fused state has too much "intrinsic inertia", imparted by 
original velocities ±P, which pushes the pulse to expand. 

At still smaller values of P, the collision also gives rise to merger of the two solitons into 
a single pulse. However, in that case, the decrease of the above-mentioned "intrinsic inertia" 
allows the fused pulse to form a stable soliton, see a typical example in Fig. [3l 
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FIG. 1: The diagram of outcomes of head-on colhsions between identical in-phase sohtons, moving 
at velocities ±P (the values of P are shown on the logarithmic scale). The black, light gray, and 
dark gray colors mark parameter regions where the merger, transition to a chaotic delocalized 
state, and quasi-elastic passage have been observed. 

Before proceeding to the presentation of results obtained for collisions between solitons 
with opposite signs, it is relevant to mention that we have also considered collisions of 
solitons with the phase difference of 7r/2. In that case (not shown here in detail), the merger 
of two solitons into a single one is also observed at small velocities, but under more specific 
conditions. In particular, for P < 0.78, the merger takes place in the region of 1 < /? < 3, 
which is essentially narrower than the merger region for in-phase solitons, cf. Fig. [H The 
reduction of the merger region is quite natural, as the interaction between the solitons is 
weaker in this case than between in-phase solitons. 
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FIG. 2: (Color online) A typical example of the transition to delocalization along the x direction, 
triggered by the merger of colliding in-phase solitons with velocities P = ±0.78, for (5 = 1. The 
evolution of the wave field is illustrated by a set of snapshots of the distribution of | U (X, T) \ . Note 
that the profiles displayed in the first panel (for Z = 0) pertain not to the established solitons, but 
to input ([2]). Profiles of the established solitons generated from these inputs can be seen below in 
the first panels of Figs. [3]and[5l 

C. Head-on collisions between solitons with opposite signs 

The diagram summarizing outcomes of collision between solitons with opposite signs 
is displayed in Fig. HI The out-of-phase solitons with large velocities ±P pass through 
each other quasi-elastically. Moreover, because the simulations were run with the periodic 
boundary conditions, we could actually observe multiple collisions, which kept their quasi- 
elastic character indefinitely long, similar to what was observed in the case of collisions 
between fast in-phase solitons, as mentioned above. At smaller (intermediate) velocities, 
collisions between the solitons with opposite signs give rise to the formation of chaotic 
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FIG. 3: (Color online) An example of the merger of colliding in-phase solitons into a single stable 
one, for P = ±0.078, /3 = 1. 

states indefinitely expanding along X, again similar to what was reported above for the case 
of in-phase solitons (therefore, examples of these outcomes are not displayed here). 

A new outcome is produced by collisions of slowly moving out-of-phase solitons. They 
do not merge into a single pulse, because this is prevented by the mutual repulsion. As seen 
in the left-hand part of Fig. [5], the solitons come close to each other and then bounce back, 
but do not separate. Instead, they arrange themselves into a persistent bound state in the 
form of a "wobbling dipole" . This observation is illustrated by the right-hand side of Fig. 
m which displays trajectories of centers of both solitons in the plane of {X,Z). Persistent 
oscillations of the bound solitons lasted as long as the simulations were run. Note that the 
range of the evolution distance in Fig. ^ Z = 15 ■ 10^, is extremely large in comparison 
with the soliton's diffraction, dispersion, and filtering lengths, which can be estimated as 
Zdiflfr ~ ^disp ~ Zfiit ~ 1 in the present situation, following the usual definitions, ^difir ~ 
^disp ~ W^/P, Zfiit ~ W^, where Wx and Wt are the soliton's widths in the spatial and 
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FIG. 4: The diagram of head-on colhsions between identical sohtons with opposite signs. As in 
Fig. m the hght and dark gray colors designate regions of the transition to chaotic delocalized 
states, and quasi-elastic passage, respectively. In the white area, the slowly colliding solitons form 
horizontal or vertical wobbling dipoles, see the text. 

temporal directions [19] (here, the notation is implied to be as in Eq. ([1])). Dependences of 
the amplitude and frequency of the oscillations on GVD coefficient (3 is not shown here, as 
the dependence is quite weak. 

It is relevant to mention that similar stable wobbling bound states of dissipative solitons 
were reported in simulations of a system of two CQ CGLEs coupled by cubic terms, in the 



ID setting, with a group-velocity mismatch between the two equations [20| (in those works, 
they were called "zigzag" states). On the other hand, in the framework of the single CQ 
CGLE in ID, truly stable bound states do not exist, although some of them may be almost 
stable 13], with the difference that the phase shift between the two solitons is 7r/2 (rather 
than vr). Thus, the existence of the stable oscillatory dipolar bound states in the single 
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FIG. 5: (Color online) Panels in the left-hand half illustrate the formation of a horizontal "wobbling 
dipole" as a result of the head-on collision between slowly moving out-of-phase solitons, for P = 
±0.078, P = 1. The panel in the right-hand half shows trajectories of centers of the two solitons, 
in the same case. 

equation is a new feature of the 2D setting. 

The formation of the wobbling dipole follows the above scenario in the region of (3 < 3. At 
larger values of the GVD coefficient (in particular, at /? = 4 and 5), the repulsive interaction 
of slowly moving out-of-phase solitons leads to their shift in the transverse direction (to 
positive and negative values of T). Eventually, they form a vertically aligned ("stacked") 
wobbling dipole, as illustrated by a set of snapshots in the left-hand side of Fig. [61 The 
vertical dipole keeps a constant vertical (i.e., temporal) separation between the two solitons 
stacked in it, which simultaneously perform persistent oscillations in the horizontal direction 
(along X), as seen in the right-hand side of Fig. [61 

D. Collisions at a finite aiming distance 

Collisions between slowly moving in-phase solitons, with a finite offset, AT, between 
their initial trajectories in the vertical direction (alias the aiming distance), also result in 
the merger into a single quiescent pulse, provided that both AT and velocities are small 
enough. However, the transition to the delocalization, similar to that shown in Fig. [21 was 
not observed at finite AT. In the same case, collisions between solitons with opposite signs 
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FIG. 6: (Color online) The same as in Fig. [5l but in the case when the collision of slowly moving 
solitons with opposite signs ends up with the formation of the vertical "wobbling dipole", at /3 = 5. 

result in a simple dynamical effect, a rebound in the vertical direction. Namely, due to the 
repulsion between the out-of-phase solitons, the value of AT increases after the collision. 
These outcome are not shown here, as they are quite obvious. 

Varying collision velocities ±P, it is possible to find a critical value of the offset, such 
that the interaction becomes negligible if AT exceeds the critical offset. For both cases of 
the in-phase and vr-out-of phase soliton pairs, the critical value is shown, as a function of 
the velocity, in Fig. [7l 

III. THREE-SOLITONS COLLISIONS 

Once the character of the two-soliton collisions has been understood, the next natural 
step is to analyze collisions between three solitons. To this end, we took triplets of identical 
solitons, with initial velocities +P, 0, —P. The systematic analysis was restricted to the 
most interesting case of the collisions between slow solitons, with P < 0.078. 

First, we consider symmetric configurations, with solitons' signs (+, +, +) and (+, — , +). 
In the former case, the outcome of the collision in the entire region where the simulations 
were run, < /? < 7, is the merger of the triplet of in-phase solitons into a single pulse, see 
an example in Fig. [HI In the latter case, the triplet of solitons with alternating signs always 
features the transition to delocalization. A noteworthy peculiarity observed in the case of 
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FIG. 7: (Color online) The critical value of the offset between initial trajectories of two solitons 
(AT in the text), in the case of non-head-on collisions. If the offset exceeds the critical value, 
the solitons effectively cease to interact. This value is shown versus the initial velocities, ±P, 
separately for in-phase and out-of-phase soliton pairs. 



(+,—,+) is formation of a transient "wobbling tri-pole" configuration, that qualitatively 
resembles horizontal wobbling dipoles generated by collisions of two solitons with opposite 
signs, see Fig. |3 Nevertheless, the "tri-pole" eventually collapses, initiating the transition 
into an expanding chaotic delocalized state. 

The collision-induced transformation of the asymmetric triplet formed by three identical 
solitons, of type (+,+,—), was investigated too. As shown in Fig. [9|, in this case, two 
in-phase solitons merge into a single one, at the first stage of the evolution (the middle 
panel in Fig. Eventually, the intermediate pair of the two pulses also demonstrates a 
transformation into a single residual soliton. 



13 



Z=0 Z=1334 




-10 -5 5 

X 



FIG. 8: (Color online) The merger of three in-phase solitons into a single one, at (3 = 1. Initially, 
edge solitons move at velocities ±0.078. 

IV. CONCLUSION 

We have undertaken the systematic analysis of collisions between two and three solitons 
in the 2D CGLE (complex Ginzburg-Landau equation) with the CQ (cubic-quintic) non- 
linearity, which may be considered as a model of large-area laser cavities, with the solitons 
representing spatiotemporal "light bullets" in it. This model, which includes the diffraction 
in the spatial direction (X) and both the GVD (group-velocity dispersion) and spectral fil- 
tering in the temporal direction, is Galilean invariant in the former direction, which makes 
it possible to create moving solitons and collide them. Outcomes of the collisions were sys- 
tematically studied by varying the GVD coefficient, f3, collision velocity, and relative sign of 
the solitons. 

In the case of collisions between two in-phase solitons, three outcomes have been iden- 
tified: the merger into a single standing soliton, transition into an expanding delocalized 
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FIG. 9: (Color online) The same as in Fig. [HI but starting with the asymmetric triplet of solitons, 
configured as (+,+,—). 

chaotic state, and quasi-elastic passage. Collisions between solitons with opposite signs at 
small velocities lead, instead of the merger, to a novel outcome - the formation of robust 
"wobbling dipoles" , which exist in two modifications: horizontal at smaller values of /?, and 
vertical ones at larger (3. In the ID setting, oscillatory bound states of dissipative solitons 
of the "zigzag" type were found in a system of two nonlinearly coupled CQ CGLEs |20|, 
but truly stable bound states in the single equation of this type do not exist. The stable 
"wobbling dipoles", especially their "vertical" variety (the one shown in Fig. [6]), is a feature 
specific to the 2D model. We have also investigated two-soliton collisions with a finite aiming 
distance, and identified its critical size beyond which the solitons cease to interact. 

Collisions between three slowly moving in-phase solitons lead to their merger into a 
single one, while three solitons with alternating signs form a transient "tri-pole" state, 
which eventually collapses into a delocalized chaotic state. Three solitons which form an 
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asymmetric configuration, of type (+,+,—), also merge into a single pulse. 
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